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Giant Meterwave Radio Telescope observations of an 
M2.8 flare: insights into the initiation of a flare-coronal 
mass ejection event 
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M. R. Kundu, S. M.White, V. I. Garaimov § 



Abstract. We present the first observations of a solar flare with the GMRT. An 
M2.8 flare observed at 1060 MHz with the GMRT on Nov 17 2001 was associated 
with a prominence eruption observed at 17 GHz by the Nobeyama radioheliograph 
^L(< and the initiation of a fast partial halo CME observed with the LASCO C2 corono- 

graph. Towards the start of the eruption, we find evidence for reconnection above 
the prominence. Subsequently, we find evidence for rapid growth of a vertical current 
' sheet below the erupting arcade, which is accompanied by the flare and prominence 

^SJ ' eruption. 
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^ . 1. Introduction 

cn 

, We have recently started observing the Sun with the Giant Meterwave 

1^ I Radio Telescope (GMRT) near Pune, India. We report here the ob- 

servation of an interesting flare-coronal mass ejection (CME) event on 
O I Nov 17 2001. The event was analyzed using data from the GMRT, 

the Nobeyama Radioheliograph (NoRH) and the Hiraiso Spectrograph 
g3 ■ (HiRAS) in Japan, the Radio Solar Telescope Network (RSTN) in Aus- 

tralia, the Mount Wilson observatory in the USA and the Michelson 
Doppler Imager (MDI) and Large Angle Spectroscopic Coronograph 
(LASCO) onboard the Solar and Heliospheric Observatory (SOHO). 

We give a brief description of the GMRT and the solar observing 
procedure in §2, followed by a multiwavelength overview of the flare- 
CME event in §3. We interpret the data in §4 and draw conclusions in 
^5. 
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2. The Giant Metrewave Radio Telescope (GMRT) 

The GMRT is located about 80km north of the city of Pune in Maha- 
rashtra, western India (latitude 19" 06 N, longitude 74° 03 E, altitude 
650 meters). It consists of thirty 45m diameter antennas spread over 
25 km. Half of these are in a compact, randomly distributed array 
of about 1 km and the rest are spread out in an approximate 'Y' 
configuration. The shortest baseline is 100m and the longest is 26 km. 
GMRT currently operates in the frequency bands around 150, 235, 
325, 610 and 1000-1450 MHz. It is among the most sensitive telescopes 
in the world in these bands and also provides seconds of arc angular 
resolution (Ananthakrishnan & Rao, 2002). Further details about the 
GMRT can be found at the URL|http://www.gmrt.ncra.tifr.res.in 

2.1. Observing the Sun with the GMRT 

These observations represent the first attempts at imaging the Sun with 
the GMRT after most of the antennas had started functioning. Our 
observing frequencies ranged from 1060 to 1076 MHz. There were 28 
antennas working through the duration of the observations presented 
in this paper. We did not track the Sun continuously during these 
observations. Instead, we repositioned the antennas at the Sun center 
every 30 minutes. Within each 30 minute interval, the phase center 
was allowed to drift. The data were then corrected for solar rotation 
during subsequent data analysis. We typically observed the Sun for 
around 20 minutes and the phase calibrator (3C298) for 6-7 minutes in 
each 30 minute cycle. Solar attenuators of 30 dB were inserted while 
observing the Sun and they were removed while recording data on the 
calibrator. Since the GMRT antennas are not yet equipped with extra 
high noise calibrators (which can generate noise corresponding to ~ 

1 sfu) we cannot place much confidence in the absolute values of flux 
observed at this time. We therefore normalized the fluxes we observed 
to the peak flux observed by the Nobeyama Radio Polarimeter (NoRP) 
at 1 GHz (http://solar.nro.nao.ac.jp/norp/html/daily/). Since the au- 
tomatic gain control (AGC) circuits were on throughout the duration 
of the observations, the flux at the flnal peak of the flare is somewhat 
supressed relative to the preceding smaller peaks (see figure 8). We used 
an integration time of 16 seconds; this is therefore the minimum time 
cadence with which we can obtain snapshot images for this observation. 
In principle, it is possible to go down to an integration time as low as 

2 seconds but no lower, because the time constant of the AGC circuits 
is of the order of 1 second. The data were processed using the standard 
NRAO/AIPS software. Snapshot images, with a dynamic range of ~10, 
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of the flaring region were obtained every 16 seconds using a restoring 
beam size of 34 X 24 in right ascension and decHnation respectively. 



3. An Overview of the event 

A long duration M2.8 flare was observed on 17 Nov. 2001 at 1060 MHz 
with the GMRT. The X-ray flare which originated in NOAA active 
region 9704 (S18E41), started at 04:49 UT, showed a smooth rise, 



GOES-10 (M2.8 FU\RE: Nov. 17 2001) 




Time (UT Hrs.) 

Figure 1. Soft X-ray lightcurves from the GOES satellites. Upper curve: 1 - 8 A. 
Lower curve: 0.5 - 4 A. The X-ray flare started at 04:49 UT, showed a smooth rise 
and peaked at 05:25 UT and exhibited a long decay period before ending at 06:11 
UT. 

peaked at 05:25 UT and exhibited a long decay period before ending 
at 06:11 UT. The X-ray fluxes observed by GOES 10 in the 1-8 A 
and 0.5-4 A bands are shown in Figure 1. Figure 2 shows the MDI 
longitudinal magnetogram of AR 9704 (left panel) and the Mount Wil- 
son sunspot drawing for this region (right panel). An examination of 
the longitudinal magnetogram shows that it is clearly more complex 
than a simple bipole. The central black polarity spot is sandwiched by 
white polarities to its northwest and southeast, suggesting the existence 
of at least two sets of loop systems. This active region was classified 
as a /3 configuration by the Solar Geophysical Data Reports on Nov 
17 2001, but its complexity was upgraded to (3^5 on Nov 19 2003, 
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Figure 2. Left Panel: The MDI longitudinal magnetogram of AR 9704. The magnetic 
complexity of AR 9704 is clearly evident. Right Panel: The Mount Wilson sunspot 
drawing for this region. The Mount Wilson sunspot drawing on Nov 17 2003 assigns 
a /3 7 (5 category to this active region. 

when it had rotated closer to disk center. The Mount Wilson sunspot 
drawing on Nov 17 2003 assigns a category to this active region. 
Notwithstanding the disparity between the classifications of AR 9704 
on Nov 17 2003, the magnetogram and the sunspot drawing provide 
sufficient evidence for the magnetic complexity of AR 9704. The 1060 
MHz GMRT images, described later, (see figure 9) show evidence for 
two sets of loop systems in this region, providing additional evidence 
for its complexity. 

The flare was accompanied by a fast halo (~ 1500 km/s) CME that 
first appeared in the C2 field of view of LASCO at around 05:30 UT 
(http://cdaw.gsfc.nasa.gov/CME_list). Figure 3 shows the lightcurve 
of the 17 and 34 GHz emission from the NoRH and a representative 
17 GHz image from the NoRH that shows the erupting prominence. 
The bright on-disk emission at 17 GHz evident in Fig. 3 is centered 
at around S18E40. The 17 GHz emission shows a characteristically 
gradual rise and decay, indicating that it is thermal in nature at these 
frequencies. On the other hand, the impulsive nature of the 200, 410, 
500 and 1060 MHz lightcurves (see figures 6-8) is evidence of their 
nonthermal nature. The 10 minute cadence movie of the 17 GHz data 
for this event available on the Nobeyama webpage (see figure 3) suggests 
that the prominence started lifting off sometime between 04:45 UT and 
04:55 UT. A detailed analysis of the prominence eruption associated 
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Figure 3. Upper panel: The upper curve is the 17 GHz hghtcurve and the lower curve 
is the 34 GHz lightcurve from the NoRH. The vertical axis gives the flux in SFU. 
Both the hghtcurves show a gradual rise and decay, implying that the radiation at 
these frequencies is thermal in nature. This radiation arises primarily from the active 
region situated at S18E40, that is evident in the lower panel. Lower panel: Image of 
the Sun at 17 GHz at 04:55 on Nov 17 2001 from the NoRH showing the erupting 
prominence on the east limb. The prominence, which was situated towards the north- 
east of the thermal active region emission evident at S18E40, erupted at ~ 04:48 UT. 
A movie of this event is available at |http:/ /solar.nro.nao. ac.jp/norh/html/10min/] 



with this event by Kundu et al. (2003) estabUshes that the rapid phase 
of the prominence hftoff started at around 04:48 UT. 

Figure 4 shows dynamic spectrum of the event from the Hiraiso 
spectrograph. The spectrum shows drifting continuum emission start- 
ing at ~ 04:35 UT at ~ 200 MHz and drifting downwards in frequency. 
The lower envelope of this drifting continuum emission represents the 
upper end of a feature that is drifting upward in the solar corona. 
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Figure 4- Dynamic spectra from the Hiraiso spectrograph. The start of the contin- 
uum emission at ~ 200 MHz at ~ 04:35 UT is clearly evident. The lower envelope of 
the drifting continuum represents the upper end of a feature that is drifting upwards. 
The height-time profile of this lower envelope is shown in figure 5. 

The HiRAS dynamic spectrum did not have data below 50 MHz for 
this event. We were therefore able to follow the lower envelope of the 
drifting continuum only until ~ 04:50 UT. An approximate curve was 
fitted to the lower envelope of the drifting continuum. Using the for- 
mulae given by Aschwanden & Benz (1995), the frequency of emission 
was related to its height in the solar atmosphere assuming that the 
observed radiation is emitted at the fundamental plasma frequency. 
Figure 5 shows the results of these computations. It shows the upward 
motion of the feature represented by the lower envelope of the drifting 
continuum (evident in figure 4) through the solar corona. 

In contrast to the slow rise and fall in the 17 and 34 GHz light 
curves, the light curves at 200, 410, 500 and 1060 MHz show an im- 
pulsive nature. In fact, the 200 MHz emission is the first sign of non 
thermal activity associated with this event. Figure 6 shows fixed fre- 
quency lightcurves at 200 MHz (lower panel) and 500 MHz (upper 
panel) from the Hiraiso spectrograph while Figure 7 shows the fixed 
frequency light curve at 410 MHz from the Radio Solar Telescope 
Network (RSTN) telescope at Learmonth, Australia. The emission at 
200 Mhz starts rising at around 04:22 UT while that at 500 MHz 
starts rising at around 04:50 UT. Figure 8 shows the 1060 MHz light 
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Figure 5. Left panel: The height-time profile of the lower envelope of the drifting 
continuum evident in the dynamic spectrum of figure 4. It represents the upper end 
of a vertical current sheet that forms beneath the erupting prominence. Right panel: 
The velocity-time profile of the same feature. The abrupt increase in velocity at 
around 04:48 UT can be inferred from the increased curvature around this time of 
the lower envelope of the drifting continuum (figure 4). 



curve derived from the GMRT observations. This emission has three 
distinct peaks. The largest peak of the GMRT 1060 MHz hghtcurve is 
normahzed to the largest peak of the 1 GHz lightcurve from the NoRP 
( http:/ / solar.nro.nao.ac.jp/norp/html/daily/ I . The general shape of the 
GMRT 1060 MHz lightcurve is very similar to the NoRP 1 GHz lightcurve. 
However, the third and largest peak in the GMRT 1060 MHz data is 
suppressed relative to the first two peaks, since the AGC circuits were 
on during the observation. As a result, the ratios peakl/peakS and 
peak2/peak3 from the GMRT 1060 MHz lightcurve are around 50% 
larger than the corresponding ratios obtained from the NoRP 1 GHz 
lightcurve. Figure 9 shows 16 s snapshots of the 1060 MHz emission 
observed with the GMRT. The images are centered at S15E46, with 
a field of view of 500" x 500" . The first peak of the 1060 MHz GMRT 
lightcurve near 04:52 UT is associated with the brightening of a banana- 
shaped NS oriented structure near the center of the field of view of the 
images. The second peak near 05:18 UT is also associated with the 
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Figure 6. Fixed frequency lightcurves at 200 and 500 MHz from the Hiraiso spec- 
trograph. The emission at 200 Mhz starts rising at around 04:22 UT while that at 
500 MHz starts rising at around 04:50 UT. 



brightening of the same structure. The third and strongest peak, at 
around 05:32 UT, corresponds to a brightening of a NE-SW oriented 
source to the south of the 'banana'. The brightness temperature of 
the radiation at each of these three peaks is < lO^K, indicating its 
nonthermal nature. Figure 10 shows 16 s snapshot images near the 
rising edge of the first peak. The images clearly show that the first 
peak of the 1060 MHz emission occured at around 04:48 UT. 

The 1060 MHz GMRT images reveal the existence of at least two sets 
of loop systems (N-S and NE-SW). Together with the magnetograms 
and sunspot drawing (see figure 2), this provides concrete evidence for 
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Figure 7. Fixed frequency light curve at 410 MHz from the Radio Solar Telescope 
Network (RSTN) telescope at Learmonth, Australia. The vertical axis is calibrated 
in SFUs above the baseline. 



the complexity of the magnetic fields in the region from which the flare 
and CME were initiated. 



4. Initiation of the fiare-CME event: analysis 

Table 1 shows the time sequence of events leading to, and after 
the CME-flare initiation, with emphasis on the nonthermal emission. 
The third and fourth columns show the estimated height in Rq from 
which the corresponding phenomenon is believed to be emanating. We 
employ the model of Aschwanden & Benz (1995) in calculating these 
heights. The third column shows the height assuming that the radio 
emission emanates at the fundamental plasma frequency (s=l) while 
the height in the fourth column is calculated assuming that the radio 
emission emanates at the first harmonic of the plasma frequency (s=2) . 
A comparison of different models for estimating the height in the corona 
from which such radio emission emanates shows that the uncertainties 
in the height estimates can range from 5 to 10% (Mann et al. 1999). The 
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Figure 8. The 1060 MHz lightcurve from GMRT observations. The three peaks are 
clearly evident. The first onset of nonthermal emission at this frequency occurs 
at 04:48 UT. The peak of this lightcurve is normahzed with respect to the peak of 
the 1 GHz Ughtcurve from the NoRP fhttp:/ /solar. nro.nao.ac.jp/norp/html/daily/ ). 
Figure 6 shows 16 s snapshot images near the three peaks of this lightcurve. 



difference between the heiglits quoted in the third and fourth columns 
is relevant only for the nonthermal radio emission at 200, 500 and 1060 
MHz. 

Masuda et al. (1994) estimate that the hard X-ray emission at the 
top of the flaring loop originates at a height of around 18,000 km. 
We have used this number as a rough estimate of the height at which 
the hard X-ray flare (lower curve, figure 1) originates. The fifth column 
gives the start time of the phenomenon, as evident from the correspond- 
ing lightcurve. We surmise that the nonthermal emission is a signature 
of magnetic reconnection, which provides sites for particle acceleration 
and consequent nonthermal emission. The start of nonthermal emission 
at a particular frequency therefore signals the start of reconnection 
processes at the height corresponding to that frequency. 

The first sign of nonthermal activity associated with this event is 
the 200 MHz nonthermal emission quoted in item A (table 1). Detailed 
measurements of the trajectory of the prominence eruption associated 
with this event made by Kundu et al. (2003) using 17 GHz data from 
the NoRH reveal that the prominence is at ~ 1.07 at 04:22 UT. 
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Figure 9. Contours of 16 s snapshots of the GMRT 1060 MHz emission. The center 
of these images is at S15E46 and the field of view is 500 x500 . The restoring beam 
is 34" X 24". The contour levels are 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9 and 1.0 
times the flux at the largest peak, which occurs at 05:33 UT. Wo show four sots 
of 16 s snapshots starting at 04:52:23 UT, 05:20:28 UT and 05:33:32 UT. These 
correspond to the throe peaks of the 1060 MHz lightcurve (figure 8). The emission 
at the peaks of 04:52 UT and 05:20 UT are from a NS oriented loop system, whereas 
the emission at the peaJi at 05:33 UT is from a NE-SW oriented loop. The brightness 
temperature of the emission at each of these peaks is < 10* K. 



The 200 MHz nonthermal emission peak of item A (table 1) occurs at 
1.08 (s=l) or 1.14 Rq (s=2), which is clearly above the prominence. 
This could be consistent with the predictions of the breakout model for 
solar eruptions (Antiochos et al. 1999), in which the first signature 
of eruption is the null point reconnection above the central arcade 
which is going to erupt. According to this model, the central arcade 
rises slowly during the initial buildup stage, due to shearing at its 
footpoints, pushing the overlying constraining field lines and distorting 
the X-type null point (above the central arcade). Beyond a certain 
point, the current layer at the overlying null point becomes sufficiently 
thin so that ideal MHD evolution is no longer possible and reconnection 
sets in, so that the overlying constraining field lines connect to adjoining 
arcades. The measurements of Kundu et al. (2003) however reveal that 
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Figure 10. Snapshots at 16 second intervals near the rising edge of the first peak in 
the GMRT light curve. The contour levels and restoring beam are the same as those 
in figure 9. It is clearly evident that the first peaJj of the 1060 MHz emission occurs 
between 04:47:32 and 04:47:48 UT. 



the prominence had started rising weh before 04:22 UT, the time at 
which the 200 MHz nonthermal emission peak of item A (table 1) 
occurs. This could be because the evolution of the null point above 
the central arcade prior to 04:22 UT was ideal, since the current layer 
was not thin enough for reconnection to set in. However, this does 
remain an open question. As mentioned earlier, it should also be borne 
in mind that there could be an uncertainty of the order of 5-10% in 
the estimate of the height from which the 200 MHz emission originates. 
The uncertainty in the projected position of the filament in the 17 GHz 
images used by Kundu et al. (2003) is much smaller; it is of the order 
of 1 beam, which is around 15 arcseconds or 0.01 R0. 

The next noteworthy feature evident from Table 1 is the near- 
simultaneity of items B, C, D, E, F and G. The 200, 410, 500 and 1060 
MHz lightcurves exhibit peaks between ~ 04:48 - 04:50 UT. The X-ray 
flare and the prominence eruption also take place in this time interval. 
Measurements of the prominence trajectory by Kundu et al. (2003) 
using NoRH 17 GHz data show that items B-G (Table 1) occur well 
below the prominence. The near-simultaneous nonthermal emission at 
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Table I. Sequence of events in CME-flare initiation 



ITEM 


DESCRIPTION 




Rq 

(s=l) 


R© 

(s=2) 


START 
(UT) 


A 

A 


zUU iVlrlZ 




1 OQ 

i.Uo 


11/1 
i.l4 


U4:zz 


"D 
D 


Lower envelope of 












drifting continuum^ 




l.iob 


i.oo 


04:48 


C 


410 MHz ^ 




1.044 


1.08 


04:48 


D 


1060 MHz N-S source 




1.02 


1.04 


04:48 


E 


Prominence eruption® 








~ 04:48 


F 


X-ray flare® 




1.02 


1.02 


04:49 


G 


500 MHz (2)'' 




1.037 


1.067 


04:50 


H 


1060 MHz N-S source 


(2)« 


1.02 


1.04 


05:13 


I 


500 MHz (2)'' 




1.04 


1.06 


05:13 


J 


1060 MHz NE-SW source^^ 


1.02 


1.04 


05:30 


K 


CME" 




5 


5 


05:30 



Notes 

^ Peak in fixed frequency lightcurve from Hiraiso, figure 6 
^ Rapid increase in velocity of lower envelope of drifting 

continuum (Figs. 4 & 5 ) 
^ Peak in fixed frequency lightcurve from Leaxmonth, figure 7 
* GMRT 1060 MHz data; figures 8 and 10 

NoRH 17 GHz data; figure 3 and Kundu et al. (2003) 
® GOES data; figure 1, lower curve, and Masuda et al. (1994) 

Peak in fixed frequency lightcurve from Hiraiso; figure 6 
® GMRT 1060 MHz data; figures 8 and 9 
^ Fixed frequency lightcurve from Hiraiso; figure 6 
° GMRT 1060 MHz data; figures 8 and 9 
^ LASCO C2 data 



heights ranging from 1.02 to 1.18 R0 could thus arise from a vertical 
current sheet below the erupting arcade. During this time interval, 
the lower envelope of the drifting continuum in the HiRAS dynamic 
spectrum (see figures 4 and 5) could represent the top of this vertical 
current sheet. The velocity of this feature increases abruptly at around 
04:48 UT (right panel of figure 5) , as is evident from the increase in the 
curvature of the lower envelope of the drifting continuum around 04:48 
UT (figure 4) . This could be taken to mean that the top of the vertical 
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current sheet accelerates upward around this time. Our observations 
thus indicate that, during the interval 04:48 - 04:50 UT: 

1. a vertical current sheet spanning at least 1.02 to 1.18 Rq exists 
(items B-G, table 1) 

2. the top of the current sheet accelerates upward (figure 5, right 
panel) 

3. the prominence eruption and X-ray flare occur 

The formation of a vertical current sheet under the erupting arcade is 
a feature of several flux rope models for CME initiation (e.g., Forbes 
1991; Lin & Forbes 2000; Mikic & Linker 1994; Amari et al. 2000). 
Reconnection at the vertical current sheet allows the eruption to pro- 
ceed in these flux rope models. The height of the flux rope undergoes a 
sudden jump as a result of the (equally sudden) formation of a current 
sheet underneath it (Lin et al. 1998; Forbes 2000). The formation of a 
vertical current sheet underneath the erupting arcade is also a feature 
of the breakout model, although it is a secondary process not integral 
to the breakout per se (Antiochos et al. 1999; Klimchuk 2000). In the 
breakout model, reconnection in the vertical current sheet results in an 
increase in the speed of the eruption. The strong heating associated 
with this current sheet reconnection is manifested as the X-ray flare. 

Items H, 1 and J (table 1) are post-eruptive phenomena. In partic- 
ular, the 1060 MHz emission from the NE-SW oriented GMRT source 
(item J) is most likely caused by flare accelerated electrons trapped in 
the NE-SW oriented loop system. This is manifested as type IV-like 
decimetric continuum emission on the Hiraiso dynamic spectrum. 

5. Conclusions 

We have presented some of the first GMRT observations of a solar erup- 
tion. We have used GMRT observations at 1060 MHz in conjunction 
with data from the NoRH, the Hiraiso solar observatory and LAS CO to 
investigate the onset mechanism of a flarc-CME event on Nov 17 2001. 
We find evidence for reconnection above the erupting arcade that could 
be construed as evidence in favor of the magnetic breakout model. We 
find evidence for the formation of a vertical current sheet below the 
erupting prominence. An abrupt increase in the velocity of the top of 
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the vertical current sheet is acompanied by the M2.8 X-ray flare and 
the prominence eruption. The features associated with the formation of 
the vertical current sheet are consistent with the predictions of several 
variants of flux rope CME initiation models, as well as the magnetic 
breakout model. 
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